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\ (57) Abstract: Disclosed is a method 

for assigning a pulse response to one of 
a plurality of pulse types with different 
decay times, particularly the pulse 
response of a scintillation detector 
with at least two scintillation materials 
having different decaying features. 
The innovative method comprises the 
following steps: - an output pulse response 
is recorded and the pulse response is 
transformed into an electric signal, the 
amplitude time response of said electric 
signal representing the pulse response 
of the output pulse; - the amplitude time 
response is transformed into the frequency 
space so as to receive an amplitude 
frequency response representing the 
output pulse; - the amplitude frequency 
response is scaled so as to receive a scaled 
amplitude frequency response; - the scaled 
amplitude frequency response is compared 
to a predefined reference response; and 
- the output pulse response is assigned to 
one of the pulse types based on the result 
of the comparison. 
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Beschreibung: 



ygrfahren zur Zuordrmna eine * PnlalaufB zu einem von ei- 
n^r Mehrzphl von Pu.i st.ypen ve rs c hiedene r fthkli.ngzeit und 
Vorrich-Knna zud sssen Durnhf uhruna 



Die Erfindung betrifft ein verfahren zur Zuordnung eines 
Pulsverlaufes zu einem von einer Mehrzahl von Pulstypen 
verschiedener Abklingzeit, insbesondere eines Pulsverlau- 
fes eines Szintillationsdetektors mit zumindest zwei 
Szintillationsmaterialien unterschiedlicher Abklingcha- 
rakteristik. Des weiteren betrifft die Erfindung eine 
Vorrichtung zur Durchfvihrung des Verfahrens. 

in der Nuklearmedizin werden radioaktiv markierte Sub- 
stanzen zur in-vivo Messung physiologischer Prozesse ver- 
wendet. Eine empf indliche nuklearmedizinische Methode ist 
die Positronen-Emissions-Tomographie (PET). Dabei werden 
isotope verwendet, die unter Aussendung eines Positrons 
zerf alien. Ein besonderer Reiz dieses Verfahrens besteht 
darin, dafl einige der physiologisch interessantesten Ele- 
mente, wie Kohlenstoff, Stickstoff oder Sauerstoff kurz- 
lebige radioaktive Isotope haben, die beim Zerfall ein 
Positron emittieren. Das ausgesandte Positron annihiliert 
im Gewebe innerhalb einer Strecke von 1-2 mm mit einem 
Elektron und erzeugt zwei, in entgegengesetzte Richtungen 
ausgesandte Gammaquanten einer Energie von 511 keV. Die 
gleichzeitige Detektion der Gammaquanten in gegenviberlie- 
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genden Detektoren, die ublicherweise in einem Detektor- 
ring angeordnet sind, bildet die Grundlage fur eine tomo- 
graphische Bildrekonstruktion. 

jeder Detektor des Rings enthalt einen Szintillationskri- 
stall, der beim Auftreffen eines Gammaquants einen Licht- 
blitz mit einer fur das Szintillationsmaterial charakte- 
ristischen Abklingzeit aussendet. Beispielsweise weist 
das oft verwendete Wismutgermanat (BSO) eine Abklingzeit 
von 300 ns auf, was zusammen mit der relativ geringen 
Lichtausbeute die erreichbare Ortsauf losung beschrankt. 
Neuere Szintillationsmaterialien mit grofiem Potential fur 
PET-Anwendungen umfassen Lutetium-Oxyorthosilikat (LSO) 
mit einer 8-fach hoheren Lichtausbeute als BSO und einer 
Abklingzeit von 40 ns, und LuAP (LuAl0 3 ) mit einer Ab- 
klingzeit von 17 ns. Die Lichtpulse der Szintillatorkri- 
stalle werden von einem Photodetektor , typischerweise ei- 
nem Photovervielfacher oder einer Photodiode erfafit, der 
zur weiterverarbeitung geeignete elektrische Signale aus- 
gibt . 

Moderne Detektoren enthalten zur Steigerung der Ortsauf- 
16sung zwei oder mehr verschiedene Szintillationskristal- 
le mit unterschiedlichen Abklingzeiten, die mit demselben 
Photodetektor gekoppelt sind. Durch eine Analyse der Aus- 
gabesignale des Photodetektors kann aufgrund der ver- 
schiedenen Abklingzeiten bestimmt werden, welcher Detek- 
torkristall angesprochen hat, und somit die Wechselwir- 
kungstiefe (depth of interaction, DOI) im Detektor erfaJJt 
werden . 
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Urn diese Zuordnung der Photodetektor-Ausgabesignale zu 
einem der Szintillationskristalle durchzuf uhren, wird ub- 
licherweise ein integrierender Verstarker beim Auslesen 
der Photodetektorsignale eingesetzt. Ein Vergleich der 
Anstiegszeit (definiert als die Zeit von 10% bis 90% des 
Maximalwerts) erlaubt dann eine Unterscheidung der ver- 
schiedenen Detektorpulse. 

In einem anderen Ansatz wird fur jedes Detektorsignal ein 
Startwert (Start des Pulses) und ein Stopwert (Ende des 
Pulses, bzw erster Nulldurchgang) erzeugt. Ein Zeit-zu- 
Amplituden-Wandler liefert dann ein Ausgangssignal mit 
einer Amplitude, die der Zeitdauer zwischen dem Start- 
und Stopwert proportional ist, und das zur Diskriminie- 
rung der Detektorpulse verwendet werden kann. 

Allen genannten Unterscheidungs- und Zuordnungsverf ahren 
ist gemein, daJ3 sie eine niedrige Zuverlassigkeit bei 
verrauschten Signalen aufweisen und eine hohe Abtastrate 
fur die zuzuordnenden Signalpulse erfordern. 

Der Erfindung, liegt daher die Aufgabe zugrunde, ein Ver- 
f ahren der eingangs genannten Art anzugeben, das eine im 
Hinblick auf Abtastrate effiziente und zuverlSssige Zu- 
ordnung eines Pulsverlaufs zu einem von einer Mehrzahl 
von Pulstypen verschiedener Abklingzeit ermoglicht. 

Diese Aufgabe wird erf indungsgemaB durch ein verfahren 
zur Zuordnung eines Pulsverlauf es zu einem von einer 
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Mehrzahl von Pulstypen verschiedener Abklingzeit, insbe- 
sondere eines Pulsverlauf es eines Szintillationsdetektors 
mit zumindest zwei Szintillationsmaterialien unterschied- 
licher Abklingcharakteristik gelost, das die Verfahren- 
schritte umfafit: 

Aufnehmen eines Ausgangspulsverlauf es und Wandeln des 
Pulsverlauf es in ein elektrisches Signal , dessen 
Amplituden-Zeit-Verlauf den Pulsverlauf des Ausgangs- 
pulses reprasentiert; 

Transf ormieren des Amplituden-Zeit-Verlauf s in den 
Frequenzraum urn einen den Ausgangspuls reprasentie- 
renden Amplituden-Frequenz-Verlauf zu erhalten; 

Normieren des Amplituden-Frequenz-Verlauf s um einen 
normierten Amplituden-Frequenz-Verlauf zu erhalten; 

Vergleichen des normierten Amplituden-Frequenz- 
Verlauf s mit einem vorbestimmten Ref erenzverlauf ; und 

Zuordnen des Ausgangspulsverlauf s zu einem der 
Pulstypen auf Grundlage des Ergebnisses des Ver- 
gleichs. 

Die Erfindung beruht damit auf dem Gedanken, die Zuord- 
nung zu einem der Pulstypen nicht im Zeitraum f sondern im 
Frequenzraum durchzufiihren. Dazu werden die gemessenen 
Pulsverlauf e nach ihrer Transformation normiert und mit 
einem Ref erenzverlauf verglichen, der gleichsam einen 
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Grenzverlauf zwischen zwei verschiedenen Pulstypen dar- 
stellt. Aus dem Ergebnis dieses Vergleichs kann dann der 
gemessene Pulsverlauf einem der Pulstypen zugeordnet wer- 
den. Dabei erfordert ein Unterscheidung zwischen zwei 
Pulstypen in der Regel einen Vergleich. Wird eine Zuord- 
nung zu einem von mehr als zwei Pulstypen vorgenommen, 
konnen auch mehrerer Vergleiche mit mehr als einem Refe- 
renzverlauf erf order lich sein. 

Die Ausfuhrung des Vergleichs im Frequenzraum kann bei 
hoherer Zuverlassigkeit wesentlich effektiver (geringere 
Abtastrate) erfolgen als im Zeitraum. Zugleich stehen ef- 
fiziente Algorithmen zur Verfugung urn den Zeitverlauf des 
gemessenen Pulses in den Frequenzraum zu trans formier en. 

In einer bevorzugten Ausgestaltung des Verfahrens wird 
das elektrische Signal vor der Transformation in den Fre- 
quenzraum einer Analog/Digital-Wandlung mit einer vorbe- 
stimmten Abtastrate unterzogen, urn einen den Ausgangspuls 
reprasentierenden , diskreten Amplituden-Zeit-Ver lauf zu 
erhalten. 

Zweckmaflig wird bei dem erf indungsgemafien Verfahren bei 
dem Schritt des Trans formierens der diskrete Amplituden- 
Zeit-Ver lauf einer diskreten Fouriertransf ormation (DFT) 
unterzogen, urn einen diskreten Amplituden-Frequenz- 
Verlauf zu erhalten. insbesondere kann die diskrete Fou- 
riertransf ormation mit einer schnellen Fouriertransf orma- 
tion (FFT) berechnet werden. 
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Bei dem Schritt des Normierens wird vorteilhaf terweise 
der Amplituden-Frequenz-Verlauf auf die Amplitude bei der 
Frequenz Null bezogen, insbesondere wird der Amplituden- 
Frequenz-Verlauf zweckmafiig durch die Amplitude bei der 
Frequenz Null geteilt. 

in einer bevorzugten Weiterbildung des erf indungsgemaflen 
Verfahrens wird bei dem Schritt des Vergleichens der Dif- 
ferenz-Verlauf des normierten Amplituden-Frequenz- 
Verlauf s und des Ref erenzverlauf s gebildet, der Diffe- 
renz-Verlauf fur jede Frequenz mit einem vorbestimmten 
Wichtungsfaktor multipliziert , und als Zuordnungsparame- 
ter die Summe dieser Produkte iiber alle Frequenzen be- 
rechnet . 

Dabei kann die Zuordnung zu einem der Pulstypen auf 
Grundlage des Vorzeichens und/oder des Betrags des Zuord- 
nungsparameters erfolgen. Insbesondere kann im Fall, wo 
der Ausgangspulsverlauf einem von zwei Pulstypen ver- 
schiedener Abklingzeit zugeordnet wird, die Zuordnung auf 
Grundlage nur des Vorzeichens des zuordnungsparameters 
erfolgen. Dies gestattet einen besonderes einfache Ent- 
scheidung iiber die Zuordnung eines gemessene Pulses zu 
einer der Pulstypen. 

Der Wichtungsfaktor ist in einer Ausgestaltung des Ver- 
fahrens derart gewahlt, daB er ein Maximum bei kleinen 
Frequenzen hat und zu grofien Frequenzen hin abnimmt. Dies 
tragt der Beobachtung Rechnung, dafi die niedrigeren Fre- 
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quenzen starker als hohere Frequenzen zur Diskriminierung 
der Pulstypen beitragen. 

Dabei ist es im vorstehend erlauterten Zusammenhang be- 
vorzugt, wenn die hochf requenten Komponenten des elektri- 
schen Signals vor der Transformation in den Frequenzraum 
durch ein Tief paflf ilter unterdriickt werden. Hierdurch 
konnen storende Ruckf altungen von Signalen, deren Fre- 
quenz hoher als die halbe Abtastrate liegt, bei der Digi- 
talisierung (ADC) des Signals vermieden werden. 

Das Verfahren kann in besonders vorteilhaf ter Weise da- 
durch weitergebildet sein, daB vor dem Aufnehmen zuzuord- 
nender Pulsverlaufe ein Ref erenzverlauf bestiinmt wird 
durch 

Aufnehmen von EichpulsverlSuf en mit bekannter Zuord- 
nung zu zwei Pulstypen mit verschiedener Abklingzeit 
und Wandeln jeweils des Pulsverlauf es in ein elektri- 
sches Signal , dessen Amplituden-Zeit-Verlauf den 
Pulsverlauf des Eichpulses reprasentiert; 

Transformieren jeweils des Amplituden-Zeit-Verlauf s 
in den Frequenzraum urn einen den Eichpuls reprasen- 
tierenden Amplituden-Frequenz-Verlauf zu erhalten; 

Normieren des Amplituden-Frequenz-Verlauf s um einen 
normierten Amplituden-Frequenz-Verlauf zu erhalten; 
und 
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Festlegen eines Ref erenzwertes fur jede Frequenz der- 
art, dafl die Amplitudenwerte der Eichpulse des ersten 
Pulstyps fur diese Frequenz im wesentlichen groBer 
als der Ref erenzwert , und die Amplitudenwerte der 
Eichpulse des zweiten Pulstyps fur diese Frequenz im 
wesentlichen kleiner als der Ref erenzwert sind. 

insbesondere kann der Ref erenzwert fur jede Frequenz mit 
Vorteil dadurch festgelegt werden, dafl jeweils ein Ampli- 
tuden-Histogramm fur die verschiedenen Pulstypen aufge- 
stellt wird, der Schnittpunkt der Einhullenden der Histo- 
gramme bestimmt wird, und der bestimmte Schnittpunkt als 
der Ref erenzwert festgelegt wird. 

Als besonders vorteilhaft wird es im vorstehend erlSuter- 
ten zusammenhang angesehen, wenn vor dem Aufnehmen zuzu- 
ordnender Pulsverlaufe der wichtungsf aktor fur jede Fre- 
quenz bestimmt wird durch 

Bestimmen des mittleren Betrags der Abweichung der 
Amplitudenwerte von dem Referenzwert fur diese Fre- 
quenz fur jeden der Pulstypen; und 

Festlegen eine hohen Wichtung fur diese Frequenz bei 
einem groBen mittleren Abweichungsbetrag und einer 
niedrigen wichtung flir diese Frequenz bei einem nied- 
rigen mittleren Abweichungsbetrag. 
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Eine erfindungsgemafle Vorrichtung zur Durchfuhrung des 
beschriebenen Verfahrens umfaBt 

ein Mittel zum Aufnehmen eines Ausgangspulsverlauf es 
und einem Mittel zum Wandeln des Pulsverlauf es in ein 
elektrisches Signal, dessen Amplituden-Zeit-Verlauf 
den Pulsverlauf des Ausgangspulses reprasentiert ; 

ein Mittel zum Trans formieren des Amplituden-Zeit- 
Verlaufs in den Frequenzraum urn einen den Ausgangs- 
puls reprasentierenden Amplituden-Frequenz-Verlauf zu 
erhalten; 

ein Mittel zum Normieren des Amplituden-Frequenz- 
Verlaufs urn einen normierten Amplituden-Frequenz- 
Verlauf zu erhalten; 

ein Mittel zum Vergleichen des normierten Amplituden- 
Frequenz-Verlauf s mit einem vorbestimmten Referenz- 
verlauf; 

ein Mittel zum Zuordnen des Ausgangspulsverlauf s zu 
einem der Pulstypen auf Grundlage des Ergebnisses des 
Vergleichs und zur Ausgabe des Zuordnungsergebnisses ; 

ein Mittel zum Bestimmen des Ref erenzverlauf es ; 
ein Mittel zum Bestimmen der Wichtungsf aktoren. 

Dabei konnen die Mittel grundsatzlich jeweils sowohl 
durch eine Hardwarekomponente, wie einen anwenderprogram- 
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mierbaren Logikschaltkreis (FPGA), durch Software in ei- 
nem geeignet eingerichteten Rechner, oder durch Firmware 
gebildet sein. Ebenso kdnnen die Mittel in einem digita- 
len signalprozessor (DSP) implementiert sein. 

weitere vorteilhafte Ausgestaltungen, Merkmale und De- 
tails der Erfindung ergeben sich aus den abhangigen An- 
spriichen, der Beschreibung des Ausfiihrungsbeispiels und 
den Zeichnungen. 

Die Erfindung soli nachfolgend anhand eines Ausfiihrungs- 
beispiels im zusammenhang mit den Zeichnungen naher er- 
lautert werden. Es sind jeweils nur die fur das Verstand- 
nis der Erfindung wesentlichen Elemente dargestellt. Da- 
bei zeigt 

Figur 1 eine Blockdiagramm einer Vorrichtung zur Zuord- 
nung eines gemessenen Pulsverlauf es nach einem 
Ausfuhrungsbeispiel der Erfindung; 

Figur 2 eine Auftragung von Amplituden-Zeit-Verlauf en 
fur je zwei Signalpulse von zwei verschiedenen 
Szintillationskristallen; 

Figur 3 eine Auftragung der Amplituden-Frequenz- 
Verlaufe der Signalpulse von Fig. 2; 



Figur 4 



eine Auftragung der normierten Amplituden- 
Frequenz-Verlaufe der Signalpulse von Fig. 2 
und eines Ref erenzverlauf s ; 
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Figur 5 eine Auftragung der Differenz zwischen den nor- 
mierten Amplituden-Frequenz-Verlauf en und dem 
Referenzverlauf von Fig. 4; 

Figur 6 ein Diagramm, das die Wichtungsf aktoren fur die 
Frequenzkomponenten von Fig. 5 angibt; 

Figur 7 eine Auftragung der gewichteten Dif f erenzver- 
laufe von Fig. 5; und 

Figur 8 eine Auftragung des nach dem erf indungsgemaBen 
Verfahren bestimmten Zuordnungsparameters fur 
100 Ausgangspulse von zwei verschiedenen Szin- 
tillationskristallen . 



Ein Beispiel einer Anwendung des erf indungsgemaBen Ver- 
fahren zur Zuordnung eines Pulsverlauf es zu zwei Pulsty- 
pen wird nachfolgend mit Bezug auf die Figuren 1 bis 8 im 
Detail beschreiben. 

Der verwendete Szintillationsdetektor 10 enthalt zwei 
verschiedene Szintillationskristall-Schichten 12 und 14, 
die hintereinander vor einem Photodetektor 16 angeordnet 
sind. Der Szintillationskristall 12 besteht im Ausfuh- 
rungsbeispiel aus LSO (Lutetium-Oxyorthosilikat ) mit ei- 
ner Abklingzeit von 40 ns, und der Szintillationskristall 
14 aus LuYAP (Lu^Y^lOa) mit einer Abklingzeit von 17 ns. 
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Zur Bestimmung des Ref erenzverlauf es und der Wichtungs- 
faktoren wurden zunachst Eichpulse mit bekannter Zuord- 
nung zu dem schnellen LuYAP-Kristall oder dem langsameren 
LSO-Kristall auf genommen . Die vom Photodetektor 16 gelie- 
ferten elektrischen Pulsverlaufe der Eichpulse wurden mit 
einem Butterworth-Filter 22 zweiter Ordnung und einer 
Grenzfrequenz von 20 MHz tief paJ3-gef iltert . Dann wurden 
die Pulse mit einer Abtastrate f s von 40 MHz, entspre- 
chend dem doppelten der Grenzfrequenz des Tief pa!3f ilters 
von dem A-D-Wandler 2 4 digitalisiert . 

Das Ergebnis nach Tief paB-Filterung und Digitalisierung 
ist in Fig. 2 fur jeweils zwei Pulse des LSO-Kristalls 12 
(Bezugszeichen 42H, 42N) und des LuYAP-Kristalls 14 
(Bezugszeichen 44H, 44N) gezeigt. Dabei ist zur besseren 
Illustration des Verfahrens jeweils ein Puis mit hoher 
Energie (LSO: 42H, LUYAP: 44H) und ein Puis mit niedriger 
Energie (LSO: 42N, LUYAP: 44N) ausgewahlt und darge- 
stellt. 

Dann werden die Amplituden-Zeit-Verlauf e der Zuordnungs- 
vorrichtung 26 eingegeben, die in anwenderprogrammierba- 
ren Logikschaltkreise 28, 30, 32, 34 und 36 umfafit, deren 
Funktion aus der nachf olgenden Beschreibung deutlich 
wird. 

Die Amplituden-Zeit-Verlaufe der Eichpulse wurden zu- 
nachst mittels einer 32-Punkt-FFT fouriertransf ormiert 
(Bezugszeichen 28), wodurch bei der genannten Abtastrate 
eine AuflSsung Af von 40 MHz/32 = 1,25 MHz erzielt wird. 
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Da sich die Pulsformen der beiden Detektoren vor allem im 
Bereich niedrigerer Frequenzen unterhalb von 10 MHz un- 
terscheiden, wurden auch bei einer niedrigeren Abtastrate 
von f s = 2 0 MHz , einem Tief paflf ilter mit einer Grenzfre- 
quenz von 10 MHz und einer 16-Punkt-FFT, bei derselben 
Auflosung Af = 1,25 MHz, gute Zuordnungsergebnisse er- 
zielt werden. Wie weiter unten im Detail beschrieben, er- 
gab der Rest mit der hoheren Abtastrate eine sicherere 
Zuordnung der Pulse von 97% erfordert aber wegen der ho- 
heren Punktzahl der FFT hoheren Aufwand fur die Berech- 
nung. 

Das Ergebnis der 32-Punkt FFT fur die vier Pulsverlaufe 
von Fig. 2 ist in Fig. 3 dargestellt. Die i-te FFT- 
Komponente auf der horizontalen Achse entspricht dabei f L 
= (i-l)*l,25 MHz , die vertikale Achse gibt die Amplitude 
A(fi) der Pulsverlaufe bei der Frequenz fi an. Wegen der 
Symmetrie der 32 FFT-Komponenten sind in der Fig. 3 nur 
die ersten 16 FFT-Komponenten gezeigt. Die Amplituden- 
Frequenz-Verlauf e der vier Eichpulse sind durch Quadrate 
52H (LSO, hohe Energie), Kreise 54H (LuYAP, hohe Ener- 
gies Sterne 52N (LSO, niedrige Energie) bzw. Kreuze 54N 
(LuYAP, niedrige Energie) dargestellt. Die unterschiedli- 
che Gesamtenergie ist an den unterschiedlichen DC- 
Amplituden bei der Komponente i=l deutlich zu erkennen. 

Zur Normierung wird der Amplituden-Frequenzverlauf jeden 
Pulses durch seine erste Frequenzkomponente A(f x ) geteilt 
(Bezugszeichen 30). Die resultierenden normierten Ampli- 
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tuden-Frequenz-Verlaufe A^fi) sind dann nicht mehr von 
der Energie der Pulse abhangig. 

Fig. 4 zeigt die normierten Verlaufe fur LSO mit Bezugs- 
zeichen 62H und 62N, fiir LuYAP mit Bezugszeichen 64H und 
64N. Wie aus der Fig. 4 deutlich zu erkennen, sind die 
Signale der verschiedenen Detektoren im Frequenzraum 
durch eine Lucke voneinander getrennt, die bei den nied- 
rigeren Frequenzkomponenten f 2 bis f B besonders ausge- 
pragt ist. Ein Ref erenzverlauf 66 , der in der Mitte der 
Lucke zwischen den beiden Pulstypen verlauft, ist eben- 
falls in Fig. 4 dargestellt. Wie die Figur ebenfalls 
zeigt , kann mehr als eine Frequenzkomponente flir die Un- 
terscheidung der beiden Pulstypen verwenden werden. 

Der Ref erenzverlauf 66 wurde fiir jede Frequenz ± L dadurch 
erhalten, da/3 die Amplitudenwerte flir eine Reihe von 
Eichpulsen jeden Pulstyps in ein Amplituden-Histogramm 
eingetragen wurden. Der Schnittpunkt der Einhullenden der 
jeweiligen Amplituden-Histogramme ergibt den Wert des Re- 
f erenzverlauf s R(fj.) fur diese Frequenz f ± . 

Nun wurde die Differenz zwischen den normierten Amplitu- 
den-Frequenz-Verlaufen A N (fi) und dem Ref erenzverlauf 
R(fi) berechnet (Bezugszeichen 32). Das Ergebnis zeigt 
Fig. 5 fiir die LSO-Pulse als 72H und 72N, und fiir die Lu- 
YAP-Pulse als 74H und 74N. Die Dif f erenzwerte sind dabei 
fur die LSO-Pulse groBtenteils positiv, fiir die LuYAP- 
Pulse im wesentlichen negativ. Wie aus Fig. 5 ersicht- 
lich, lassen sich die verschiedenen Pulse bei niedrigeren 
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Frequenzkomponenten gut trennen , f lir hohere Frequenzkom- 
ponenten geht die Unterscheidungskraf t jedoch gegen Null. 

Urn diesem Umstand Rechnung zu tragen, wurde fur jede Fre- 
quenzkomponente ein Wichtungsf aktor VI (f L ) festgelegt, der 
ein Mal3 fiir die Eignung dieser Frequenzkomponente zur 
Pulstypen-Diskriminierung darstellt. Dabei wird jeweils 
der mittlere Betrag der Abweichung der Amplitudenwerte 
von dem Referenzwert flir jeden der Pulstypen als MaJ3 flir 
die Unterscheidungseignung genommen. Aus den in Fig. 5 
dargestellten Amplituden lassen sich etwa die in Fig. 6 
gezeigten Wichtungsf aktor en WffJ fur jede Frequenzkompo- 
nente f ± bilden. Wie auch anschaulich aus Fig. 5 klar, 
erhalten dabei die Komponenten mit niedrigem Index, ins- 
besondere die Frequenzkomponenten bei 2,5 MHz und 3,75 
MHz, ein hohes Gewicht, wahrend die Wichtungsf aktoren flir 
die nicht zur Diskriminierung beitragenden hoheren Fre- 
quenzen auf Null gesetzt sind. 

Der so ermittelte Ref erenzverlauf R(fi) und die Wich- 
tungsf aktoren W(f ± ) werden gespeichert, und konnen dann 
flir die Zuordnung gemessener Pulsverlaufe verwendet wer- 
den, deren Ursprung von einem der beiden Detektoren nicht 
bekannt ist, wie nachfolgend dargestellt. 

Zur Illustration wird die Zuordnung anhand der bisher 
verwendeten Eichpulsen des LSO und des LuYAP Kristalls 
erlautert. Wird ein neuer Puis unbekannter Zuordnung ge- 
messenen, so werden die oben beschriebenen Schritte mit 
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Ausnahme der Berechnung des Ref erenzverlauf s und der 
Wichtungsf aktoren durchgef uhrt . 

Dann wird die Differenz zwischen dem normierten Amplitu- 
den-Frequenz-Verlauf A N (fi) und dem Ref erenzverlauf R(fi) 
fiir jede Frequenz mit dem Wichtungsf aktor multipliziert 

D w (fJ = W(fJ * (A^fJ - R(f A )) 

urn eine gewichtete Differenz D w (fi) zu erhalten 
(Bezugszeichen 34). Die so gewichteten Werteverlauf e sind 
in der Fig. 7 flir die LSO-Pulse als 82H bzw. 82N und fiir 
die LuYAP-Pulse als 84H und 84N dargestellt. 

Ein Aufsummieren liber die Beitrage von alien Frequenzen 
ergibt dann den Zuordnungsparameter Z flir jeden Puis 
(Bezugszeichen 34). Zu Zuordnung gemigt das Vorzeichen 
von Z. Ein positiver Wert von Z f iihrt zu einer Einordnung 
des gemessenen Pulses als Signal des LSO-Kristalls , ein 
negativer Wert zu einer Einordnung als Signal des LuYAP- 
Kristalls (Bezugszeichen 36). 

Zur Untersuchung der Zuverlassigkeit des erf indungsgema- 
J3en Verfahrens wurden nach einraaliger Festlegung des Re- 
f erenzverlauf s R(f ± ) und der Wichtungsf aktoren W(f ± ) ein- 
hundert Ausgangspulse, fiinfzig des LSO-Kristalls und 
fiinfzig des LuYAP-Kristalls nach dem vorstehenden erlau- 
terten Verfahren zugeordnet. Figur 8 zeigt jeweils den 
Zuordnungsparameter Z fiir die LSO-Pulse (Bezugszeichen 
92) und die LuYAP-Pulse (Bezugszeichen 94). Die Z-Werte 
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fur die vier Eichpulse sind als 92H, 92N, 94H und 94N 
ebenfalls eingezeichnet . Wie aus der Figur ersichtlich, 
flihrt das erf indungsgemafie Verfahren bei 97% der Pulse zu 
einer korrekten Zuordnung. 
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Anspruche; 



Verfahren zur Zuordnuna eines Pulsverlauf es zu einem von 
einer Mehrzahl von Pulstypen verschiedener Abklinqzeit 
und Vorrichtuna zu dessen Durchfuhruna 



1. Verfahren zur Zuordnung eines Pulsverlauf es zu einem 
von einer Mehrzahl von Pulstypen verschiedener Ab- 
klingzeit, insbesondere eines Pulsverlauf es eines 
Szintillationsdetektors mit zumindest . zwei Szintil- 
lationsmaterialien unterschiedlicher Abklingcharak- 
teristik, umfassend die Verf ahrenschritte : 

- Aufnehmen eines Ausgangspulsverlauf es und Wandeln 
des Pulsverlauf es in ein elektrisches Signal , des- 
sen Amplituden-Zeit-Verlauf den Pulsverlauf des 
Ausgangspulses reprasentiert ; 

- Trans formieren des Amplituden-Zeit-Verlauf s in den 
Frequenzraum um einen den Ausgangspuls reprasen- 
tierenden Amplituden-Frequenz-Verlauf zu erhalten; 

- Normieren des Amplituden-Frequenz-Verlauf s um ei- 
nen normierten Amplituden-Frequenz-Verlauf zu er- 
halten; 

- Vergleichen des normierten Amplituden-Frequenz- 
Verlauf s mit einem vorbestimmten Ref erenzverlauf ; 
und 
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- Zuordnen des Ausgangspulsverlauf s zu einem der 
Pulstypen auf Grundlage des Ergebnisses des Ver- 
gleichs . 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet , 
daB das elektrische Signal vor der Transformation in 
den Frequenzraum einer Analog/Digital-Wandlung mit 
einer vorbestimmten Abtastrate unterzogen wird, um 
einen den Ausgangspuls reprasentierenden , diskreten 
Amplituden-Zeit-Verlauf zu erhalten. 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, 
daJ3 bei dem Schritt des Transf ormierens der diskrete 
Amplituden-Zeit-Verlauf einer diskreten Fourier- 
transformation (DFT) unterzogen wird, um einen dis- 
kreten Amplituden-Frequenz-Verlauf zu erhalten. 

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, 
daJ3 die diskrete Fouiriertrans formation mit einer 
schnellen Fouriertransf ormation (FFT) berechnet 
wird. 

5. Verfahren nach einem der vorigen Anspruche, dadurch 
gekennzeichnet, da/3 bei dem Schritt des Normierens 
der Amplituden-Frequenz-Verlauf auf die Amplitude 
bei der Frequenz Null bezogen wird, insbesondere da!3 
der Amplituden-Frequenz-Verlauf durch die Amplitude 
bei der Frequenz Null geteilt wird. 
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6. Verfahren nach einem der vorigen Anspruche, dadurch 
gekennzeichnet, dai3 bei dem Schritt des Vergleichens 
der Dif f erenz-Verlauf des normierten Amplituden- 
Frequenz-Verlauf s und des Ref erenzverlauf s gebildet 
wird, der Dif f erenz-Verlauf fur jede Frequenz mit 
einem vorbestimmten Wichtungsf aktor multipliziert 
wird, und als Zuordnungsparameter die Summe dieser 
Produkte iiber alle Frequenzen berechnet wird. 

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet , 
da/3 die Zuordnung zu einem der Pulstypen auf Grund- 
lage des Vorzeichens und/oder des Betrags des Zuord- 
nungsparameters erf olgt . 

8. Verfahren nach Anspruch 6 oder 7, dadurch gekenn- 
zeichnet , daJ3 der Ausgangspulsverlauf zu einem von 
zwei Pulstypen verschiedener Abklingzeit zugeordnet 
wird, und die Zuordnung auf Grundlage nur des Vor- 
zeichens des Zuordnungsparameters erf olgt, 

9. Verfahren nach einem der Ansprliche 6 bis 9, dadurch 
gekennzeichnet, daB der Wichtungsf aktor ein Maximum 
bei kleinen Frequenzen hat und zu groBen Frequenzen 
hin abnimmt. 

10. Verfahren nach einem der vorigen Ansprliche, dadurch 
gekennzeichnet, daB hochf requente Komponenten des 
elektrischen Signals vor der Transformation in den 
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Frequenzraum durch ein Tief pa!3f ilter unterdruckt 
werden. 

11. Verfahren nach einem der vorigen Anspruche, dadurch 
gekennzeichnet , da/3 vor dem Aufnehmen zuzuordnender 
Pulsverlaufe ein Ref erenzverlauf bestimmt wird durch 

- Aufnehmen von Eichpulsverlauf en mit bekannter Zu- 
ordnung zu zwei Pulstypen mit verschiedener Ab- 
klingzeit und Wandeln jeweils des Pulsverlauf es in 
ein elektrisches Signal, dessen Amplituden-Zeit- 
Verlauf den Pulsverlauf des Eichpulses reprasen- 
tiert; 

- Transformieren jeweils des Amplituden-Zeit- 
Verlaufs in den Frequenzraum um einen den Eichpuls 
reprasentierenden Amplituden-Frequenz-Verlauf zu 
erhalten; 

- Normieren des Amplituden-Frequenz-Verlauf s um ei- 
nen normierten Amplituden-Frequenz-Verlauf zu er- 
halten; und 

- Festlegen eines Ref erenzwertes flir jede Frequenz 
derart, dafl die Amplitudenwerte der Eichpulse des 
ersten Pulstyps fiir diese Frequenz im wesentlichen 
groJier als der Ref erenzwert, und die Amplituden- 
werte der Eichpulse des zweiten Pulstyps fiir diese 
Frequenz im wesentlichen kleiner als der Referenz- 
wert sind. 
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12. Verfahren nach Anspruch 11 , dadurch gekennzeichnet , 
daJ3 der Referenzwert fur jede Frequenz dadurch fest- 
gelegt wird, da!3 jeweils ein Amplituden-Histogramm 
fiir die verschiedenen Pulstypen aufgestellt wird, 
der Schnittpunkt der Einhullenden der Histogramme 
bestimmt wird, und der bestinunte Schnittpunkt als 
der Referenzwert festgelegt wird. 

13. Verfahren nach Anspruch 11 oder 12, dadurch gekenn- 
zeichnet, daJ3 vor dem Aufnehmen zuzuordnender Puls- 
verlaufe der Wichtungsf aktor fur jede Frequenz be- 
stimmt wird durch 

- Bestimmen des mittleren Betrags der Abweichung der 
Amplitudenwerte von dem Referenzwert fur diese 
Frequenz fur jeden der Pulstypen; und 

- Festlegen eine hohen Wichtung fiir diese Frequenz 
bei einem groflen mittleren Abweichungsbetrag und 
einer niedrigen Wichtung fiir diese Frequenz bei 
einem niedrigen mittleren Abweichungsbetrag. 

14. Vorrichtung zur Durchfiihrung des Verfahrens nach ei- 
nem der Anspruche 1 bis 13 , mit 

- einem Mittel zum Aufnehmen eines Ausgangspulsver- 
laufes und einem Mittel (10) zum Wandeln des Puls- 
verlaufes in ein elektrisches Signal, dessen 
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Amplituden-Zeit-Verlauf den Pulsverlauf des Aus- 
gangspulses reprasentiert ; 

- einem Mittel (28) zum Transf orraieren des Amplitu- 
den-Zeit-Verlauf s in den Frequenzraum um einen den 
Ausgangspuls reprasentierenden Amplituden- 
Frequenz-Verlauf zu erhalten; 

- einem Mittel (30) zum Normieren des Amplituden- 
Frequenz-Verlauf s um einen normierten Amplituden- 
Frequenz-Verlauf zu erhalten; 

- einem Mittel (32,34) zum Vergleichen des normier- 
ten Amplituden-Frequenz-Verlauf s mit einem vorbe- 
stimmten Ref erenzverlauf ; und 

- einem Mittel (36) zum Zuordnen des Ausgangspuls- 
verlaufs zu einem der Pulstypen auf Grundlage des 
Ergebnisses des Vergleichs und zur Ausgabe des Zu- 
ordnungsergebnisses . 

15. Vorrichtung nach Anspruch 14 , bei dem das Mittel 
(32 f 34) zum Vergleichen und das Mittel (36) zum Zu- 
ordnen durch einen anwenderprogrammierbaren Logik- 
schaltkreis (FPGA) oder einen digitalen Signalpro- 
zessor (DSP) oder (PC) gebildet ist. 

16. Vorrichtung nach Anspruch 15, bei dem das Mittel 
(28) zum Transf ormieren und das Mittel (30) zum Nor- 
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mieren durch einen anwenderprogrammierbaren Logik- 
schaltkreis (FPGA) oder einen digitalen Signalpro- 
zessor (DSP) oder (PC) gebildet ist. 
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